Josephson-Majorana cycle in topological single-electron hybrid transistors 
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Charge transport through a small topological superconducting island in contact with a normal and 
a superconducting electrode occurs through a cycle which involves coherent oscillations of Cooper 
pairs and tunneling in/out the normal electrode through a Majorana bound state, the Josephson- 
Majorana cycle. We illustrate this mechanism by studying the current-voltage characteristics of two- 
terminal superconductor - topological superconductor - normal metal single-electron transistors. At 
low bias and temperature the Josephson-Majorana cycle is the dominant mechanism for transport. 
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Introduction — Majorana bound states are zero-energy 
states that occur at the boundary or in the vortex core 
of topological superconductors jTHSj. Besides the fun- 
damental interest in asserting the existence of Majorana 
fermions, their discovery in condensed matter will be a 
fundamental step towards the realization of a topological 
quantum computer 0]. For these reasons an intense re- 
search has started to find physical systems that support 
Majorana excitations. A quantum wire with spin-orbit 
coupling in close proximity to a superconductor and in 
the presence of an external magnetic field is among the 
most promising proposals put forward so far [^{T]- 

The presence of Majorana bound states in topological 
superconductors has a number of distinct signatures. It is 
responsible for a fractional Josephson effect [3 El [HI [9] , it 
leads to a resonant Andreev reflection |in| and to anoma- 
lies in interference experiments , it can be detected by 
measuring the local density of states [121. 

In the presence of quantum fluctuations, affecting the 
phase coherence of the topological superconductors, new 
effects appear due to a parity constraint linking the dy- 
namics of the superconducting phase and of the Majo- 
rana fermions. A paradigmatic situation to explore these 
phenomena, in the setup of Refs. [SHZ], is when the su- 
perconducting island in contact with the wire is small 
enough so that charging effects [13] come into play [H]. 
Fu |15| first pointed out that the parity constraint leads 
to non-local effects in electron transport. The role of 
charging effects on the fractional Josephson effect and 
on the Coulomb blockade through a topological super- 
conducting transistor was studied in Refs. [T^ and [T7], 
respectively. 

In this work we show how the parity constraint leads 
to a new charge transport mechanism: The Josephson- 
Majorana cycle. It takes place when a topological su- 
perconducting island is coupled to a superconducting 
and to a normal lead (see Fig. [l]). Charges can flow 
through the island due to the combined effect of the 
coherent oscillations of Cooper pairs in the island and 
the tunneling between the Majorana state and the nor- 
mal leads. Although the process bears some similari- 
ties to the Josephson-quasiparticle mechanism present in 
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FIG. 1. The system consists of a superconducting grain cou- 
pled to a quantum wire connected to a normal metal and to a 
superconductor. Parameters are chosen so that the wire hosts 
two Majorana bound states at the two ends. 



Cooper pair transistors [18H20] there are important dif- 
ferences which we will discuss in the concluding part of 
the paper. The observation of transport mediated by the 
Josephson-Majorana cycle introduced here can be an ev- 
idence for the existence of Majorana fermions in systems 
proposed as topological superconducting wires. 

The model — The system we consider is illustrated in 
Fig. [T] It consists of a topological superconducting island 
(a nanowire in close proximity to a superconducting is- 
land) tunnel-coupled to a normal and a superconducting 
electrode. The island hosts two Majorana bound states 
at the ends of the wire associated with the Majorana 
operators -fi {i = L,R), = S^j and -fi = -fj. 

At low energies (much smaller than the superconducting 
gap) , the island couples to the superconducting lead only 
through Josephson tunneling. The coupling to the nor- 
mal lead occurs only via the Majorana bound state at 
the end of the island connected to the normal lead. The 
mesoscopic scale of the central island considered here re- 
quires to take into account the charging energy of the 
island. The Hamiltonian of the system reads: 



H^Hm+ Hch +Hj +nL+HT 



(1) 



where the five terms describe coupling of Majorana 
states, island charging energy, Josephson coupling to 
the superconducting electrode, metallic normal electrode, 
and tunneling to the normal lead, respectively. The 
coupling between the Majorana fermions is given by 
Hm ~ i^lLlR where the overlap A between them is 
exponentially small for distances much larger than the 
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superconducting coherence length and we will drop this 
term in the following. Two Majorana states can form a 
zero-energy fermionic level described by an annihilation 
operator d — (7^ + i^r) / which can be either occu- 
pied or empty. Therefore not only the number of excess 
Cooper pairs N but also the occupation of the d-level 
enters the charging energy of the island, 



Hch=Ec{2N + nd-ngy^ 



(2) 



where Ud — d}d counts the occupation of the d-level, 
Ug is the gate charge which can be varied continuously 
by changing the gate voltage, and Eq = e^/(2C) is the 
charging energy expressed in terms of the total capaci- 
tance C of the island. For later convenience we will label 
the eigenstates of the charging energy with the notation 
\N,nd)- In order to describe the Josephson coupling to 
the superconducting electrode we use the effective Hamil- 
tonian 



-Ej cos {(fis - (fi) 



(3) 



where Ej is the Josephson coupling energy and ips and 
if are the phase of the superconducting electrode and the 
island condensate, respectively. For the present purposes 
any possible coupling to the Majorana state can be safely 
ignored '21]. The normal lead Hamiltonian is given by 
the noninteracting model Hl — X]fc ^kcl.Ck with creation 
(annihilation) operator cj. (cfe) corresponding to a spin- 
less fermion in the free particle state k with energy ek 
inside the normal lead. Finally, the tunneling through 
the Majorana state takes the form [T7] . 



(4) 



where tk is the hopping amplitude to the fc-state in the 
lead. The advantage of this formulation introduced by 
Zazunov et al. |17j is the fact that it includes automat- 
ically the constraint on the Hilbert space, linking the 
occupation of the Majorana bound state to the parity of 
the superconducting condensate. A non-trivial dynamics 
in the problem arises since (p and N are canonically con- 
jugated variables, [N, e*"^] = e^'^. There are two terms in 
the tunneling Hamiltonian, corresponding to regular and 
anomalous tunneling. The first one describes the transfer 
of an electron from the d-level to the normal lead, and 
the second one the annihilation of a Cooper pair inside 
the island accompanied by the creation of two electrons, 
one in the d-level and one in the normal electrode. 

Josephson- Majorana cycle — At bias voltage and tem- 
perature smaller than the superconducting gap, quasi- 
particle tunneling is suppressed and coherent Joseph- 
son (Cooper pair) tunneling is necessary for transport 
through the transistor. From Eq. ([2| it follows that the 
resonance condition for coherent Cooper pair oscillations 
between two charge states that differ by one Cooper pair 
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FIG. 2. Josephson-Majorana cycle for rig ~ 1. The states 
differing by one Cooper pair, jO,0) and |1,0), are coupled via 
the Josephson term. The coupling between |0, 1) and |0,0) is 
given by the normal term in the tunneling Hamiltonian, the 
connection between by |0, 1) and jl,0) by the anomalous one. 



is fulfilled at integer values of the gate charge Ug. Around 
Ug ~ 1, e.g., the lowest energy states are |A^, Ud) — |0, 0), 
|0, 1) and |1,0). A Josephson-Majorana cycle involv- 
ing these three states is illustrated in Fig. [2] Starting 
from |0, 1), a regular tunneling process releases an elec- 
tron into the normal lead. Then, the island is recharged 
with an extra Cooper pair provided by the Josephson 
coupling. A second (anomalous) tunneling which annihi- 
lates a Cooper pair in order to create an electron inside 
the normal electrode and another one filling the d-level, 
completes the cycle. (A cycle with the reverse direction 
can be obtained by the conjugates of each tunneling pro- 
cess.) Since the energy of the state |0, 1) is lower than 
the ones connected by Josephson tunneling, there will be 
a threshold voltage for the onset of current. 

In the remainder of this work we support the simple 
picture sketched above with more detailed calculations. 
It is sufficient to compute the current- voltage characteris- 
tics to second order in the tunneling amplitudes tk- The 
average of the current operator, I = ie J^ki'^l'^k^ ^ 
— 2elm ^ J, ifcC^((i -I- e~'^'^d^), can be conveniently ex- 
pressed in terms of the reduced density matrix p of the 
topological superconducting island which is obtained af- 
ter tracing out the fermionic degrees of freedom of the 
normal metal. The average steady-state current is 

eF 



-Di^D*i^{2-F:i-F^i)]p^ra (5) 



where Dnm aud pnm are the matrix elements of the 
operators D = d + e^^'^d^ and p in the basis defined 
by the eigenstates \ipn) of the Hamiltonian Hch + Hj 
with eigenvalues £„. Furthermore, Fnm — /(A„m) ~ 
(i/n) Re^[l/2 + i(A„„ - eV) / (271 ksT)], where A„„, = 
£n~£m, /(e) = l/[l-|-e('-'=^)/'=«^] is the Fermi function 
of the normal electrode, and ^E* is the digamma tunc- 
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FIG. 3. Current-voltage characteristics for ksT — Ec /lOO, 
r = ksT/lO, Ej = Ec/5. 



tion. Finally F — 2Tr\T\'^N{ep) is the tunneling rate with 
t]^ ^ T assumed constant close to the Fermi energy e^, 
and N{ep) is the density of states in the normal metal. 

A convenient way to represent and compute the re- 
duced density matrix, in particular when higher-order 
tunneling processes are taken into account, is to use a 
real-time diagrammatic technique that has been devel- 
oped to describe transport through a metallic single- 
electron transistor [531 US]- The time evolution of the 
reduced density matrix p{t) can be cast in the form fM\ 



Pnm ~t- i^nynPnra ^ ^ 



dt'W, 



n7n,n'7n' 



{t,t')pn'm'it') . 



Expanding to second order in the tunneling, this yields 
for the steady-state limit {p — 0) [5S], 



i^nmPrmi — 2 ^ y [-Dn'n-^in'miFrnn' ~^ Fnrn') 
n' ,rn' 

~^-Dnn'-D„-^„^r{2 — -F„/,„ — Fm'n)] Pn'm' 



FIG. 4. Current-voltage characteristics and its dependence on 
the bias charge ny = CV/e and the gate charge rig = CVg/e. 
Parameters are: UbT = Ec/IQQ, T = knT/lQ, Ej = r/%/3. 
The threshold voltage to observe a current is ny — ±0.5. The 
slight displacement of the resonance condition from rig = 1 is 
due to the Lamb shift from the normal lead. 



For positive voltage ^ > and 1/2 < Ug < 3/2, the 
steady-state current is 



2Ej 



F2 + 4f 2 - 



(/) = eF <^ 
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(8) 



where £ = ^/\4E^ 
^ x^)/16, and x = 



j)]^ + Ej,g,ing) = {21 + rSx- 



x" — X")/ W, ana X = 4£'c(l — ng)/£. We immediately 
see that the current vanishes when the Josephson cou- 
pling goes to zero. The width of the resonance peaks is 
5ng = Ej/AEc and the peak value for V > becomes 
maximal for F — \J 1E'^ — Ej. 

The effect of energy renormalization terms (Lamb 
shifts) that were neglected in the analytical formulas pre- 
sented above, is highlighted in Fig. |4] where the resonance 
condition is slightly displaced from the value Ug — 1. We 
remark that the current is antisymmetric under the trans- 



[Dm'iF>*^iFin + £';„i(l — F*i)] Pnm' ■ (6) formation {V, Ug) — >■ {—V, 2 — Ug). In order to see this we 
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In Fig. [3] we show the current as a function of the gate 
charge Ug — CVg/e and the bias charge ny — CV/e. 
The current increases stepwise (smeared by temperature) 
with the bias voltage and shows resonance peaks at inte- 
ger values of Ug. In the range of the gate charge where 
only the three charge states |0, 0), |0, 1) and |1, 0) are in- 
volved, a simple analytical solution can be obtained for 
zero temperature when neglecting the imaginary parts of 
Fnm, which are associated with energy renormalization 
induced by the coupling to the bath. We find the two 
threshold voltages 



ey± = Er 



1±2 



^{1 - UgY + [Ej/AEcf 



(7) 



note that the transformation results in the substitutions 
Fnm ^-F.^n and Dn,n From Eqs. ^ and ^ 

we conclude that p remains unchanged while the current 
changes sign. 

Finally, we comment on the abrupt suppression of the 
current for large bias voltage \nv\ > 3/2 that is visi- 
ble in Fig. [4] Beyond the threshold ny = 3/2, the new 
charge state |1,1) becomes available from |1,0) by tun- 
neling. The island is then trapped in this new state since 
the Josephson coupling to the state |0, 1) is suppressed 
because of the large energy difference of |1, 1) and |0, 1). 
As a consequence, transport is blocked. 

As already anticipated in the introduction there are 
crucial differences between the Josephson-Majorana cy- 
cle and the Josephson-quasiparticle cycle of Cooper pair 
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transistors. The mechanism discussed in the present pa- 
per stems from the constraint hnking the occupation of 
the Majorana bound state with the parity of the super- 
conducting condensate. In the Josephson-quasiparticle 
cycle, it is a process related to the breaking of a Cooper 
pair into two quasi-particles with energy ^ A. Even the 
presence of an accidental zero energy state close to the 
normal lead would not give rise to the cycle described 
in this paper. A transport mechanism involving single 
particles and Cooper pairs would appear either in higher 
orders in tunneling or, again, above the superconducting 
gap. The presence of a zero-energy bound state, without 
the non-local character of the Majorana states, is not 
sufficient to realize the cycle introduced in the present 
paper because it would not display the anomalous tun- 
neling term in the Hamiltonian. 

Conclusions — In topological hybrid transistors charge 
transport is dominated at low voltages by a Josephson- 
Majorana cycle. This is a process in which the Rabi oscil- 
lations of Cooper pairs are accompanied by the tunneling 
of electrons from/to Majorana bound states. We showed 
the basic principles of this cycle in a superconductor - 
topological superconductor - normal metal transistor and 
discussed how the features of the current-voltage charac- 
teristics are linked to the presence of Majorana bound 
states. Because of the two time scales associated to the 
cycle, 1/Ej and l/F, interesting features are expected to 
appear in the frequency dependence of the current noise. 
Furthermore, more complex current cycles can be real- 
ized if more wires are deposited on top of the central 
superconducting island. The case of two wires is of par- 
ticular interest since it is the minimal system to realize 
a topological qubit even in the presence of interaction as 
discussed in Ref. HH 
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